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We calculate the thermoelectric power (or thermopower) of many semiconducting single wall 
carbon nanotubes (s-SWNTs) within a diameter range 0.5-1.5 nm by using the Boltzmann transport 
formalism combined with an extended tight-binding model. We find that the thermopower of s- 
SWNTs increases as the tube diameter decreases. For some s-SWNTs with diameters less than 
0.6 nm, the thermopower can reach a value larger than 2000/iV/K at room temperature, which 
is about 6 to 10 times larger than that found in commonly used thermoelectric materials. The 
large thermopower values may be attributed to the one-dimensionality of the nanotubes and to 
the presence of large band gaps of the small-diameter s-SWNTs. We derive an analytical formula 
to reproduce the numerical calculation of the thermopower and we find that the thermopower of 
a given s-SWNT is directly related with its band gap. The formula also explains the shape of 
the thermopower as a function of tube diameter, which looks similar to the shape of the so-called 
Kataura plot of the band gap dependence on tube diameter. 

PACS numbers: 79.10.-n,72.20.Pa,65.80.-g 


I. INTRODUCTION 

In recent years, there has been significant interest in re¬ 
search on thermoelectric phenomena due to the increase 
in the demand for alternative energy sources. Espe¬ 
cially, since thermoelectric phenomena could transform 
heat currents into electric power, thermoelectric power 
generators can perhaps be used to convert waste heat 
into electric energy for use in environmentally friendly 
applications [1 -3]. It is thus necessary to find a good 
thermoelectric material with a high thermoelectic energy 
conversion efficiency, characterized by the so-called ther¬ 
moelectric figure of merit, ZT = S 2 ctk,~ 1 T 1 where S is 
the Seebeck coefficient, also known as the thermoelectric 
power (thermopower), a is the electrical conductivity, k 
is the thermal conductivity, and T is the absolute tem¬ 
perature of the material. Over the past six decades it 
has been challenging to obtain ZT values exceeding 2, 
because the parameters of ZT are generally interdepen¬ 
dent [2, 3]. A theoretical study in 1993 predicted that 
the ZT value of low-dimensional structures could be sig¬ 
nificantly enhanced, thanks to the quantum confinement 
effect to create sharp features in the density-of-states 
(DOS) [4], This prediction was confirmed experimen¬ 
tally in 1996 using PbTe/Pbi-^Eu^Te, which exhibited 
a ZT value up to about five times greater than that of 
the corresponding bulk value [5]. It is thus intriguing 
to evaluate other low-dimensional structures that might 
have excellent thermoelectric performance, either theo¬ 
retically or experimentally. 

As a one-dimensional material, single wall carbon nan¬ 
otubes (SWNTs) were considered promising for thermo- 
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electric materials due to their novel electronic properties 
which depend on their geometrical structure [6-8]. How¬ 
ever, it has been difficult to obtain an ensemble of indi¬ 
vidual SWNTs with a specific (n, m) structure to reveal 
the precise knowledge of the dependence of the thermo¬ 
electric power of individual SWNTs on band gap and 
diameter. Most thermoelectric measurements were per¬ 
formed on bundled SWNT samples whose geometrical 
and electronic structures are complex [6-9], and thus the 
potential thermoelectric properties might have been lost 
because of interactions between different tubes [6]. The 
ZT values reported for bundled SWNTs have remained in 
the range of 10 -3 to 10 -4 [9, 10], in contrast to the com¬ 
mercial thermoelectric materials with ZT « 1 [11, 12]. 
Such bundled SWNT samples consist of a collection of 
SWNTs with different diameters, metallicities, and chi¬ 
ralities, parameters to which the electronic structure is 
very sensitive [13]. The small ZT value of the bundled 
SWNTs were mainly attributed to their low thermopower 
and high thermal conductivity, which might be a result 
of the mixture of different SWNTs and impurity in low 
concentration in the samples. 

In this work, we will focus on evaluating the ther¬ 
mopower theoretically for many SWNTs, especially in 
the case of semiconducting SWNTs (s-SWNTs), and thus 
to maximize the SWNT thermopower and to suggest a 
new route for obtaining a larger ZT for SWNTs. By 
calculating the thermopower of all individual s-SWNTs 
within a diameter range 0.5 < dt < 1.5 nm, we will show 
that, for tube diameters less than 0.6 nm under low dop¬ 
ing, the thermopower of s-SWNTs can be as large as 
2000 /TV /K at room temperature, which is large enough 
compared to the thermopower of bundled SWNTs, which 
is about 100 — 200^(V/K [6, 8, 14]. From this result, we 
believe that there is still much room available to improve 
the ZT of SWNT samples. For a more practical pur- 



2 



7"hot 


FIG. 1. (Color online) Schematic model of a thermoelectric 
device using two identical s-SWNTs, one with p-type and the 
other with n-type doping. The temperature gradient between 
the two edges of each nanotube generates an electric current. 


to the electron-phonon interaction. In particular, Jiang 
et al. showed that the relaxation time from the electron 
scattering with the TW phonon mode is independent of 
the electron energy [17]. Therefore, here we make the 
assumption that the thermopower from the Boltzmann 
transport equation can be obtained by applying the re¬ 
laxation time approximation (RTA) and we may even 
treat the relaxation time as a constant. Under the RTA, 
the thermopower or Seebeck coefficient S is expressed by 


S = - 


VU 

VT 


1 j v{e)T{e)v(e) 9 ^ ^ [e - n]g{e)de 
qT f v{e)t(e)v{e) g{s)d£ 


(1) 


pose, we also give an analytical formula to reproduce 
our numerical calculation of the s-SWNT thermopower, 
which forms a map of the s-SWNT thermopower. The 
calculated thermopower map could be useful for obtain¬ 
ing information on the s-SWNT chirality with a desired 
thermopower value and thus it offers promise for using 
specially prepared s-SWNT samples to guide the direc¬ 
tion of future research on the thermoelectricity. 

This paper is organized as follows. In Sec. II, we give 
the theoretical methods employed in this study to calcu¬ 
late the thermopower. In Sec. Ill, we discuss the ther¬ 
mopower obtained from the numerical calculation and 
compare it with the analytical formula. We then summa¬ 
rize the results and give the future perspective in Sec. IV. 
We also provide some appendices for a detailed derivation 
of the thermopower analytical formula. 

II. MODEL AND METHODS 

To utilize the s-SWNTs as a main material in future 
thermoelectric devices, we consider a model shown in 
Fig. 1, in which two identical s-SWNTs, one with p- 
type and the other with n-type doping, are connected 
in parallel. Each s-SWNT should maintain its electronic 
charge distribution in the nonequilibrium state, for ex¬ 
ample, by a temperature gradient along the tube axis. 
By having their temperature gradient VT from an edge 
of each s-SWNT to its other edge, charge carriers (elec¬ 
trons or holes) will flow with velocity v from the hot edge 
with temperature T^ ot to the cold edge with temperature 
Tcoid- The carrier distribution /o, which depends on the 
electronic energy e and chemical potential /i, is modified 
as a function of e, following the Boltzmann transport for¬ 
malism. Within such a process, an electric voltage W 
can be generated. It is also known from earlier studies 
that the electron-phonon interaction is the main factor 
determining the electrical conductivity of SWNTs [15- 
17], in which the so-called twisting (TW) phonon mode 
with a long wavelength gives the dominant contribution 


where q = ±e is the unit carrier charge, T = (Thot + 
T C oid)/2 is the average absolute temperature, v(e) is the 
carrier velocity, g{s) is the electronic (DOS), and t(e) is 
the carrier relaxation time. 

We employ both numerical and analytical methods to 
obtain S from Eq. (1). In the full numerical approach, we 
can use the BoltzTraP code [18], which is a widely-used 
package to calculate some thermoelectric properties, such 
as the thermopower and electrical conductivity. A neces¬ 
sary input for the BoltzTraP code is the electronic energy 
dispersion e(fc) for all bands (multiband structure). The 
BoltzTraP code also adopts a constant t, whose plausi¬ 
bility in the case of s-SWNTs has been justified above. 
While the BoltzTraP code is actually sufficient for ob¬ 
taining the thermopower from Eq. (1), we cannot discuss 
the physics of the thermopower of s-SWNTs without hav¬ 
ing an explicit formula for the thermopower that depends 
on some physical parameters, such as the SWNT band 
gap and geometrical structure. Therefore, we also solve 
Eq. (1) analytically by considering the valence band and 
the conduction band closest to the Fermi energy, known 
as the two-band model [19, 20]. The derivation of the 
analytical formula is explained in detail in Appendices 

a-d'. 

As the input for the BoltzTraP code, we calculate the 
energy dispersion s(k) within the extended-tight bind¬ 
ing (ETB) model developed in our group [21]. The ETB 
model takes into account long-range interactions, SWNT 
curvature corrections, and geometrical structure opti¬ 
mizations, which are sufficient to reproduce the exper¬ 
imentally observed energy band gaps of the SWNTs [21 
23]. The SWNT structure in our notation is denoted by 
a set of integers (n, m) which is a shorthand for the chi¬ 
ral vector Ch = n&i + m& 2 , where ai and &2 are the 
unit vectors of an unrolled graphene sheet [13]. The chi¬ 
ral vector C/j defines the circumferential direction of the 
tube, giving the diameter d t - Another vector perpen¬ 
dicular to C h defines the tube axis, which is called the 
translational vector T [13]. The chiral and translational 
vectors thus represent the tube unit cell. In the Boltz¬ 
TraP calculation, we use a 20 nm x 20 nm x |T| supercell, 
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FIG. 2. (Color online) Thermopower as a function of chem¬ 
ical potential and temperature for an (11,0) s-SWNT. 


where |T| (in nm) is the length of the translational vec¬ 
tor. A large supercell length in the x- and //-directions is 
chosen so as to guarantee the individual SWNTs are well- 
separated. Since the thermopower in the BoltzTraP code 
is expressed in terms of a tensor [18], the corresponding 
thermopower tensor component for a given s-SWNT is 
S zz , which is the thermopower along the tube axis direc¬ 
tion. Other tensor components are negligible. 


FIG. 3. (Color online) Thermopower as a function of chemi¬ 
cal potential for (11, 0), (12,4), and (15, 5) s-SWNTs at 300 K. 
Solid lines are obtained from the numerical calculation based 
on Eq. (1) while dashed lines are obtained from the analytical 
formula given in Eq. (2). 


final form of Scnt can written as 


Scnt = 


k B 

e 


k B T 


E g 

2 k B T 


3 

2 


E g /k B T + 3 \ 

e 2^/k B T + 1 J > 


III. RESULTS AND DISCUSSION 

In Fig. 2, we show a first example of the thermopower 
calculation result for an (11,0) s-SWNT. The ther¬ 
mopower (S zz ) is plotted versus chemical potential and 
temperature. We see that the thermopower is higher at 
the lower temperature because S oc 1/T in Eq. (1). The 
maximum thermopower obtained for the (11,0) SWNT is 
about 1420 /iV/K, which is already large for a purely in¬ 
dividual s-SWNT compared to that for bundled SWNTs 
with S of around 100—200 /iV/K [6, 8]. Next, we can also 
plot the thermopower at a specific temperature to see the 
chemical potential dependence of the thermopower. In 
Fig. 3, we show the thermopower versus chemical poten¬ 
tial for three different s-SWNT chiralities: (11,0), (12,4), 
and (15,5), at T = 300K. The solid lines in Fig. 3 
represent the numerical results. For all chiralities, the 
optimum value of the thermopower, indicated by a max¬ 
imum (minimum) along the negative (positive) axis of 
the chemical potential, arises due to the p-type (n-type) 
characteristics of the s-SWNTs, which is consistent with 
a recent experimental observation [8]. The dependence of 
the thermopower on the chemical potential implies that 
it is possible to tune the thermoelectric properties of s- 
SWNTs by applying a gate voltage, giving p-type and 
n-type control over the thermopower. 

To better understand the numerical results of ther¬ 
mopower, we have derived an analytical formula for the 
thermopower within the two-band model [19, 20]. We de¬ 
note this analytical formula of the thermopower as Scnt 
( see Appendices A-D for the detailed derivation). The 


where e is the elementary electric charge, k B is the Boltz¬ 
mann constant, and E g is the SWNT band gap. The 
Eg values adopted in Eq. (2) are obtained from previous 
ETB results [21]. The dashed lines in Fig. 3 represent 
the fit of the numerical results of the thermopower using 
Eq. (2) for three different s-SWNT chiralities. The an¬ 
alytical formula [Eq. (2)] fits the numerical results near 
/r = 0. In particular, the two optimum thermopower val¬ 
ues (maximum and minimum for p-type and n-type dop¬ 
ing, respectively) can be well-reproduced in that region, 
which implies that the energy bands near the Fermi en¬ 
ergy give the strongest contribution to the thermopower 
of s-SWNTs. The analytical results deviate from the nu¬ 
merical results at larger |/i| far from the optimum ther¬ 
mopower because the two-band model is no longer valid 
at a higher doping level. However, for the discussion in 
this paper, the two-band model is already sufficient to de¬ 
scribe the thermopower of s-SWNTs since we will mainly 
focus on the optimum values of the thermopower. 

For a more rigorous argument, we determine a condi¬ 
tion to obtain an optimized chemical potential y 0 pt from 
Eq. (2), which satisfies dScNT{Ho P t)/dn = 0. We then 
obtain 



where the + and — signs define the n-type and p-type 
contributions, respectively. From Eq. (3), we can say 
that the fj, opt values will move more distant from /i = 0 
as E s becomes larger than k B T, as shown in Fig. 4(a). 
However, due to the presence of the logarithmic term, 
/Apt is very slowly changing as a function of E g when E g 
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FIG. 4. (Color online) The optimized chemical potential 
Mopt plotted as a function of s-SWNT band gap. In panel (a), 
we scale the chemical potential and the band gap by ksT/2 
and fcflT, respectively, as described by Eq. (3). In the case of 
(b), we set a constant T = 300 K and vary E g , while in (c) 
we set a constant E g = 0.913 eV, which is the band gap value 
of an (11,0) s-SWNT, and vary the temperature. 


is much larger than fc^T. This behavior can be seen in 
Fig. 4(b), in which we show the E g dependence of /.i opt . 
For the dt range of 0.5-1.5nm, the s-SWNTs have E g 
values of about 1.58 eV down to 0.46 eV. In this case, 
Eg is about 17-61 times larger than fcgT for T = 300 K. 
With those E g values, we then obtain 0.046 < |/x opt | < 
0.062eV at a constant T = 300K [see Fig. 4(b)], which 


implies that the change in /r opt in this case is only about 
16 meV although the change in E g is as large as about 
1.12 eV for the same d t range. At room temperature, 
controlling the doping level or the chemical potential is 
thus useful to give us the optimum thermopower for the 
s-SWNTs under consideration. On the other hand, by 
decreasing T for a given E g , we can also decrease /r opt , 
as shown in Fig. 4(c), which reduces the doping level 
required to obtain the optimum thermopower. It should 
be noted that in Fig. 4(c) we intentionally set a constant 
Eg = 0.913 eV for simplicity although the s-SWNT band 
gaps in the realistic case may decrease as a function of 
temperature by about 3% when we increase T from 200 K 
to 800 K [24], 

Using both the numerical calculation by BoltzTraP 
and our analytical formula S'cnt, it is now possible for us 
to plot the thermopower of s-SWNTs over a broad range 
of dt by taking the optimum value of the thermopower. In 
the case of the analytical formula, we define the optimum 
thermopower from Eqs. (2) and (3), as follows 

^cnt = -S'cnt (m = Mo P t)- (4) 

In Figs. 5(a-b), we show the optimum thermopower val¬ 
ues of many s-SWNTs with 0.5 < d t < 1.5 nm compared 
with their corresponding band gaps as a function of di¬ 
ameter. In Fig. 5(a), we plot the optimum thermopower 
calculated from the BoltzTrap simulation (denoted by 
circles) and from SqjJt (denoted by plus symbols) on the 
same scale. We can see that the two methods show a 
good agreement. From Fig. 5(a), the thermopower of s- 
SWNTs is also found to increase as the tube diameter 
d t decreases. For some s-SWNTs with d t < 0.6 nm, such 
as those with 2 n + m = 13, i.e. the (5,3) and (6,1) s- 
SWNTs, the thermopower can reach a value more than 
2000 p,V/K. These thermopower values are about 6-10 
times larger than those found in common thermoelectric 
materials [12, 25-28]. 

The larger thermopower for smaller-diameter s- 
SWNTs can be explained by the relation of S'cnt with E g 
as shown in Eq. (2) and by the fact that E g oc 1/dt [29]. 
The one-dimensional character of the SWNT electronic 
DOS may also enhance the thermopower [4, 5]. Here, 
we should note that the thermopower of s-SWNTs as a 
function of diameter shows the nanotube family pattern, 
in which the different SWNTs with the same 2 n + m can 
be connected and they make a clearly distinct branch for 
mod(2n+m, 3) = 1 and mod(2 n+m, 3) = 2, known as the 
nanotube SI and SII family branches, respectively [29]. 
This behavior is very similar to that found in the band 
gap as a function of diameter shown in Fig. 5(b), which 
is often referred to as the Kataura plot [23, 29, 30]. This 
result also suggests that the measurement of the ther¬ 
mopower of a single chirality s-SWNT sample might be 
able to predict an exact band gap value of the s-SWNT. 
In fact, the band gap is directly connected to the ther¬ 
mopower as can be seen in the S'cnt formula [Eq. (2)]. 

Finally, we would like to briefly discuss the issues of 
maximizing the thermoelectric power factor, which is the 
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FIG. 5. (Color online) (a) Optimum thermopower val¬ 
ues for all s-SWNTs within the diameter range of 0.5-1.5nm 
plotted as a function of SWNT diameter. The temperature is 
set constant at 300 K. Numerical results from BoltzTraP are 
denoted by circles, while analytical results from Eqs. (2)-(4) 
are denoted by plus symbols, (b) The Kataura plot show¬ 
ing the family pattern of the SWNT band gap as a function 
of diameter. Solid lines are a guide for the eyes, connecting 
SWNTs with the same family number 2n + m. The SI and SII 
tubes correspond to the SWNTs having mod(2n + m, 3) = 1 
and 2, respectively. 


numerator term in the ZT formula. There are two main 
issues to which we have to pay attention. First, we may 
argue that, for s-SWNTs as a thermoelectric material, it 
might still be impossible to obtain a large ZT or a us¬ 
able device at a low doping level despite the fact that the 
optimum thermopower values are obtained near /z = 0. 
The reason is that the electrical conductivity er can be 
very small near /i = 0. This fact is also reflected in the 
conductivity equation as a function of /z [see Appendix B, 
Eq. (B4)]. However, compared to the bulk materials, the 
one-dimensional materials such as s-SWNTs have smaller 
effective mass m *, which may enhance the electrical con¬ 
ductivity due to the relation of er oc (m*) -1 / 2 , as can also 


be seen in Eq. (B4). Second, we may worry that, as we go 
to smaller diameter s-SWNTs (in which the thermopower 
is optimized), the electrical conductivity will instead be 
too small to maximize the power factor. However, we 
note that there is also a chirality dependence which could 
enhance the electrical conductivity through the effective 
mass relation. As mentioned before, a smaller to* will 
give a larger er, and thus s-SWNTs which have both small 
diameters and small to* might be useful as a thermoelec¬ 
tric material even at relatively low doping levels. 


IV. SUMMARY 

We have shown the theoretically predicted behavior 
of the thermopower of many s-SWNTs within a diam¬ 
eter range of 0.5-1.5nm. We derive a simple formula 
to calculate the thermopower of s-SWNTs from their 
band gap, which enables us to predict the optimum ther¬ 
mopower values. The optimum thermopower value of 
an individual s-SWNT (p-type or n-type) can be larger 
than 2000 /zV/K at room temperature for diameters less 
than 0.6 nm, such as the (5,3) and (6,1) s-SWNT. Our 
results highlight potential properties of small diameter 
s-SWNTs as a one-dimensional thermoelectric material 
with a giant thermopower. With the recent advances 
in the fabrication methods for specific small diameter s- 
SWNTs [14, 31, 32], we expect that the further potential 
development of s-SWNT thermoelectric devices could be 
realized in the near future. 
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Appendix A: Thermopower of nondegenerate 
semiconductors 

Here we derive a general formula for the thermopower 
of nondegenerate semiconductors as a starting point be¬ 
fore deriving the analytical formula of S'cnt [Eq. (2)]. In 
the calculation of the thermopower, we assume that the 
single wall carbon nanotubes (s-SWNTs) are nondegen¬ 
erate semiconductors. The thermopower or the Seebeck 
coefficient for a nondegenerate semiconductor can be cal¬ 
culated by solving the Boltzmann transport equation un¬ 
der the relaxation time approximation, which leads to the 
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following expression [ 20 ]: 


1 J v&T^vie) 9 -^ [e - n]g(e)de 
qT J v(£)t{e)v(£) d ^ £ ' ) g(e)d£ 


(Al) 


where q 1 e, T, and /i are the unit carrier charge, elec¬ 
tronic band energy, temperature, and chemical potential, 
respectively. The variables u(e), r(e), /o(e), and g{£) are 
the band carrier velocity, carrier relaxation (scattering) 
time, Fermi-Dirac distribution function, and the density 
of states (DOS) per unit volume, respectively, defined by 


v 2 {e) = 


2 e 

m*d' 


(A2) 


t{e) = T 0 £ r , (A3) 


fo(e) 


1 

1 _|_ e {s-fi)/k B T ’ 


(A4) 


9(e) 


1 (^L\ d/ \d/2-i 

2 7 3-d 2 d-l 7r d/2 r (|) \ H 2 J 


(AS) 


where d = 1, 2,3 denotes the dimension of the material, 
to* is the effective mass of electrons or holes, r is a char¬ 
acteristic exponent, tq is the relaxation time constant, 
and L is the confinement length for a particular material 
dimension. Substituting Eqs. (A2)-(A5) into Eq. (Al) 
yields 


( 

V 



(A 6 ) 


To simplify Eq. (A 6 ), we define the following variables: 
the reduced band energy £ = e/(fcsT), the reduced chem¬ 
ical potential 77 = /i/(/csT), and the Fermi-Dirac inte¬ 
gral Fj(rj) = J £ J '/o(£)d£- Inserting these quantities into 
Eq. (A 6 ) gives 


S = (n - ! +r + 1 x Jd/2+r_\ 

9 ( f + r F d / 2 +r-i J 


(A7) 


Since (£ — 77 ) > 3 for nondegenerate semiconductors, we 
can use an approximation of Fj(rj) ss + 1), where 

r(j) is the gamma function, to obtain 


ks 

9 


( l+r + 1 T(l+r + l) \ 

{ l+r nl+r) J' 


(AS) 


Using the recursion formula T(j + 1) = jT(j), the ther¬ 
mopower of nondegenerate semiconductors within the 
one-band model can be written as 


S 


ks 

9 


77 - r — 1 

2 


This last equation is still insufficient to derive the ther¬ 
mopower of s-SWNTs since the s-SWNTs are consid¬ 
ered as nondegenerate semiconductors with two energy 
bands. In this case, we also need an expression of electri¬ 
cal conductivity because the semiconductor within the 
two-band model includes a conduction band for elec¬ 
trons and a valence band for holes following the formula 
S = ( a n S n + a p S p )/(a n + a p ) [19], where S niP and <t„, p 
are, respectively, the thermopower and electrical conduc¬ 
tivity for the n-type or p-type semiconductors. The ex¬ 
pression specifying the electrical conductivity for a single 
energy band is derived in Appendix B. 


Appendix B: Electrical conductivity for 
nondegenerate semiconductors 

The electrical conductivity is expressed as [20] 

cr = -q 2 J v(£)t(£) v(e) g(e)de. (Bl) 

Substituting Eqs. (A2)-(A5) and the Fermi-Dirac inte¬ 
grals into Eq. (Bl) yields 

_ 2 < 7 2 t 0 (f + r) 1 /2 m*\ d/2 

m*d L3-d 2 <i-i 7r d/2 r (|) J 

x ( k B T) d ' 2+r Fd/ 2+r-l- (B2) 

By applying the approximation Fj(rj) ss e n T'(j -)- 1) for 
nondegenerate semiconductors, we can write the electri¬ 
cal conductivity, 

2<? 2 t 0 (| + r) _1_ /2m*\ d/2 

° ~ m*d L 3 ~ d 2 d ' 1 7r d / 2 r(|) \ h 2 ) 

x (fc B T) d / 2 +VT(J+r^ , (B3) 

which finally becomes 

fj = Vrod+r) (fc B T) d / 2 +W(f+r) 

d L 3 - d 2 d l 2 ~ 3 Tx d l 2 h d Y{\) K ’ 

(B4) 

We will use Eq. (B4) for calculating the electrical con¬ 
ductivity to derive the thernropower of two-band semi¬ 
conductors in the next section. 


Appendix C: Thermopower of two-band 
semiconductors 

The thermopower of two-band semiconductors is de¬ 
fined by [19] 


S = 


U n. + (J P S P 


(A9) 


(Cl) 
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tion of the SWNT diameter. SI and SII tubes correspond to 
the SWNTs having mod(2n + m, 3) = 1 and 2, respectively. 
Solid lines connect SWNTs with the same 2 n + m value. 


where cr„ jP and S n>p are expressed as 

2 q 2 r 0 (l+r) ( k B T) d / 2+r T{± + r) 
an,P ~ d L3-d 2 rf/2- 1 7r d/2^r(|) 

x (C2) 

and 

S n , P = (v n, P - - r - 1^ , (C3) 


respectively. Substituting Eqs. (C2) and (C3) into 
Eq. (Cl), and after doing some algebra, we can obtain 



<Jp 

(C4) 

where a n /a p = (m*/m*) d ^ 2 1 e r,n ~ rtp = Ae rin ~ Vp , with 

A = (m^/m*)^ 2 1 . Here we have r] n — r) p = 2r] p 
and r] n + ?y p = -r] g , where = n/(k B T) and rj g = 
Eg/(k B T). The thernropower of nondegenerate semicon¬ 
ductors within the two-band approximation can then be 
written in terms of i) p , rj g , r, d , and A as 


S = 


k B 

e 



Tj g -f- 2 v -\- d + 2 \ 
Ae 2r i p + 1 ) ' 


(C5) 



FIG. 7. (Color online) The percentage error, or the discrep¬ 
ancy between the analytical and the numerical results of the 
thermopower calculations for each s-SWNT, is plotted versus 
the SWNT diameter. The discrepancy increases linearly with 
increasing the SWNT diameter, as indicated by the fitted 
dashed lines. 


effective masses m* in the s-SWNTs can be calculated 
using the effective mass formula m* = h 2 (d 2 e/dk 2 )~ 1 , 
where e(k ) is the electronic energy dispersion within the 
extended tight binding (ETB) model [21]. We can obtain 
A as a function of diameter, as can be seen in Fig. 6, in 
which we show A within a diameter range of 0.5-1.5 nm. 
In this diameter range, we have A « 1. With such an 
approximation, and also assuming that the carrier relax¬ 
ation time is constant [which gives r = 0 according to 
Eq. (A3)], the thernropower of s-SWNTs is then given 
by 


S'CNT 


^b_(„ _ Vs _ 3 , Vg + 3 \ 

e V ^ 2 2 + e 2r > p + 1 ) ' 


(Dl) 


The thernropower can be rewritten in terms of /r and E g 
as 


<Scnt — 


k B 

e 


k B T 


E g 3 E g /k B T - 1"3\ 

2 k B T ~ 2 + e 2 ^/ k B T + i J ' 

(D2) 


where k B is the Boltzmann constant and E g is taken from 
the ETB calculation [21]. We see that Eq. (D2) is nothing 
but Eq. (2). In this derivation, the reason why we put 
r = 0 is that the electron relaxation time r in s-SWNTs 
is determined mainly by the electron-phonon interaction 
with the TW phonon mode (see the main text, Sec. II), 
where the relaxation time is taken to be independent of 
the electron energy [17]. Therefore, we can write r = r 0 
or equivalently r = 0. 


Appendix D: Thermopower of s-SWNTs 

We now finally have all the information needed to de¬ 
rive Scnt- Since s-SWNTs are one-dimensional, we have 
d = 1 and A = (m*/m*) . The electron and hole 


Appendix E: Comparison between numerical and 
analytical methods 

To verify the accuracy of the Sq([ t in fitting the nu¬ 
merical results of the s-SWNT thermopower, we show 













in Fig. 7 the difference of the thermopower obtained 
from the analytical and numerical calculations in terms 
of the error percentage. This error percentage variable 
is the difference in the thermopower calculated by using 
the 5° p n % formula with respect to the numerical results 
for each s-SWNT diameter. We obtain the error values 
ranging from —2% to 4% for both p-type and n-type s- 


SWNTs. The error values increase with the increase of 
the tube diameter because E g oc 1/dt and also because 
the formula for S'cnt [Eq. (2)] was derived by assuming 
s-SWNTs as nondegenerate semiconductors. Therefore, 
larger band gaps or smaller diameter s-SWNTs should be 
more accurately fitted by our S'cnt approximation. 
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